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We develop a theory for the field emission effect in Lut- 
tinger liquids and single-wall carbon nanotubes at the level 
of the energy resolved current distribution. We generalise 
Fowler-Nordheim relations. Just below the Fermi edge, we 
find a power-law vanishing current distribution with the den- 
sity of states exponent. The current distribution above the 
Fermi edge owes its existence to a peculiar interplay of in- 
teractions and correlated tunnelling. It displays a non-trivial 
power-law divergence just above the Fermi energy. 

PACS numbers: 03.65.X, 7I.10.P, 73.63. F 



The field emission (FE) effect has attracted consid- 
erable attention for a long time. Energy resolved cur- 
rent densities, otherwise also termed total energy distri- 
bution functions (TEDs), have been thoroughly studied, 
both experimentally and theoretically (for a review see 
0). A theoretical framework for the FE was devised in 
the late 1920's by Fowler and Nordhcim |]. They de- 
scribe the emission process as electron tunnelling from a 
solid into vacuum through an asymmetrical triangularly- 
shaped barrier. The height of the barrier with respect to 
the Fermi level is given by the work function W , and the 
slope on the vacuum side is proportional to the electric 
field F applied to the electrode, as shown in Fig. [l[ 




FIG. 1. Schematic representation of the field emission 
process. Here D is the conduction band-width of the emit- 
ter, W is the work function and V is the bias voltage for a 
related tunnelling problem (see text). In the inset: Diagram 
describing the scattering of a hot hole Pi into Pf and creation 
of an electron-hole pair e — pc with subsequent tunnelling of 
the electron into the right electrode (dashed line). 



If the tunnelling amplitude is small, the energy re- 



solved current J{ijj) is proportional to the probability 
n(w) for an electron to have energy w, to the energy de- 
pendent transmission coefficient I'(w), and to a factor T 
responsible for the tip geometry: 



J(w) = TV{uj)n{uj) . 



(1) 



In the noninteracting case n(w) is given by the Fermi 
distribution function. The (quasi-classical) transmis- 
sion probability for a triangular barrier is given by [^i 
2?(w) ~ exp (-4V2^(VF - ujf^/ifiF), where m is the 
electron mass and the electron energy ld is from now on 
measured relative to the Fermi energy Ep- At zero tem- 
perature the emerging spectrum has a sharp threshold 
at the Fermi edge (no particles with energies above the 
edge) and is essentially constant in the vicinity of the 
edge. Integrating over all energies one obtains the one 
dimensional version of the Fowler-Nordheim formula re- 
lating the full current to the electric field's strength [^,0. 

What is even more exciting, measurements by Lea and 
Gomer revealed that in real systems there is a finite 
current above the Fermi edge ('secondary current'). They 
found a singular TED approximately of the form 



(2) 



This can only be explained by taking into account the 
interplay between the electron-electron interactions and 
correlated tunnelling as there is only one way for the 
'secondary' electron to gain additional energy. That is 
by scattering on hot holes left by the 'primary' electron 
after tunnelling into the vacuum. Gadzuk and Plum- 
mer applied Boltzmann-type equations for the particle 
balance to explain this phenomenon Thereby they 
used the scattering cross-section of the process depicted 
in Fig. |l| and obtained the TED consistent with (||). 

In this letter we investigate the FE effect in Luttinger 
liquids (LLs). The LL model is the basic model for the 
description of one-dimensional (ID) interacting electrons 
It is relevant for the physics of organic quasi-lD sys- 
tems and quantum wires, for recent reviews see |7|j^. Our 
main interest is due to the fact that the LL model is also 
responsible for low-energy properties of single-wall car- 
bon nanotubes (SWNT) as discussed in a series of papers 
Q over the last few years. Especially the conductance 
properties of SWNTs in various setups were exhaustively 
investigated both theoretically and experimentally [ p^ . 
Apart from recent widespread use of nanotubes as STM 
tips, they are expected to be applied in electronics as field 
emitters in high-resolution displays and cathode tubes 
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pl| . While there are several experimental investigations 
of the FE from nanotubes the theory of this pro- 
cess in LLs was barely discussed in the literature so far. 
In this letter we undertake to close this gap, as far as 
SWNTs are concerned. From the theoretical point of 
view the most interesting part of the theory is that of 
the secondary current. To our knowledge the latter has 
not yet been studied experimentally in SWNTs. So, af- 
ter a brief but comprehensive discussion of the primary 
current from LLs we turn to the TED above the Fermi 
edge. 

In order to apply Eq. (|^) to LLs one has to know the 
TED n(u;) in the emitter. In computing the primary 
current one neglects higher-order tunnelling processes, so 
that n{uj) is a Fourier component of the equilibrium cor- 
relation function (V'^(i)V'(O))- Here 'ip'^{t) is a local elec- 
tron creation operator at the tip of the electrode x — 0. 
The TED in this case is proportional to the local den- 
sity of states (LDOS) at the tip of the emitter. Using the 
standard bosonisation scheme for open-boundary LLs , 
we write 

^{x = 0,t) = (2^ao)"'/' exp[iq^{x - 0, t)/^ , (3) 

where ao is a lattice constant and g is the LL parameter. 
The gaussian chiral bosonic field 0(a;, t) is governed by 
the LL Hamiltonian 

and is periodic in 2L, where L is the system size. (From 
now on we set h — vp = e = 1.) For a half- infinite 
system the LDOS's energy dependence is known [0,D, so 
that the TED can readily be written as 



n{uj) = e{-uj)\uj\^/s-i/aoD^/sY(l/g) 



(4) 



where F stands for the gamma function and D is the 
band-width of the LL. Plugging this into Eq. (||) one 
observes two important facts: (i) at the lowest order in 
tunnelling the TED above the Fermi energy is still zero, 
even for the interacting system, and (ii) below the Fermi 
energy the TED has a power-law singularity. Integrating 
over all energies we establish the Fowler-Nordheim-like 
formula for LLs: 



J 



1/29 



exp 



4fcf ' 



(5) 



Now we generalise these results for SWNTs. These sys- 
tems are known to be described as four-channel LLs . 
Three channels (^c- (charge-flavour), (total spin), 
(/)<,- (spin-flavour) are non-interacting. The fourth chan- 
nel <i)c+ (total charge, or the plasmon mode) possesses 
the LL parameter K ~ (\ + AUq/tt)^'^/'^ , where Uq is the 
zero Fourier component of the screened Coulomb poten- 
tial. Both the experimental pOl] and the theoretical 0] 



estimates place K somewhere in between 0.15 and 0.3. 
Note that though we now have four channels the field- 
operator actually factorises as [§[ 



^ ~ exp{i0c+/(4V/s:) + i(0c- + (t>s+ + 0s-)/4} , (6) 

Therefore also the correlation functions factorise and the 
above results, (||) and (^, are valid for SWNTs given the 
substitution 
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(7) 



In the rest of this paper we investigate the secondary 
current above the Fermi edge. The first challenge we face 
is to put the theory of on modern footing. To do so, 
let us take a broader view of the problem and consider 
tunnelling between two conductors at a finite bias volt- 
age V, see Fig. |l|. Clearly this becomes equivalent to the 
FE problem when the bias is large, or when the energy uj 
of the omitted electron is small as compared to V . (From 
now on we concentrate on the vicinity of the Fermi edge 
and neglect all trivial energy dependences of the tun- 
nelling probabilities.) Allowing for finite bias is not only 
beneficial technically (as it makes some integrals con- 
verge) but has important physical consequences. Namely, 
qualitative results for the secondary current will also be 
valid for usual non-equilibrium tunnelling contacts be- 
tween conductors (more about it elsewhere |13 ). The 



model Hamiltonian for the problem can thus be taken as 

i7 = i/o[c]+-ffLLM+7[^^(0)c(0)+ct(0)7A(0)] , (8) 

where {x) is the electron creation operator on the right 
of the barrier and 7 is the effective tunnelling ampli- 
tude (7^ ~ 'D(Ep)). (In our approach, it is impor- 
tant for the electron system on the receiving end to be 
non-interacting, at least asymptotically at large distances 
from the barrier.) For simplicity we use a ID model with 
a single emission channel, this will not affect the expo- 
nents. 

Since we are dealing with a profoundly non-equilibrium 
phenomenon we employ the Keldysh technique . The 
Keldysh Green's functions are given by 

g{x,x';t-t') = ~i{Tcijix,t)^\x\t')), (9) 
G(x, x']t - t') = -i{Tcc{x, t)c\x', t')) , 

where Tq stands for the Keldysh contour ordering oper- 
ator Q. The TED is related to the following Keldysh 
component 



n{Q,uj) = -ig-+(0,0;w) 



(10) 



The zero-order perturbative expansion of the TED in 
the tunnelling amplitude 7 results in already discussed 
Eq. (^). The next non- vanishing terms are of order 
7^ and, as pointed out in Refs. contain the sec- 

ondary emission processes. These secondary electrons 
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result from inelastic scattering of hot holes, left by the 
primary electrons, see Fig. |l|, and possess energies higher 
than the Fermi energy psf . There are various diagram- 
matic contributions to Eq. (^0|) which describe such pro- 
cesses. However, after a closer look at the scattering 
process shown in Fig. |l| one realises that in the relevant 
diagrams the created particles are annihilated in the same 
way they were created. The outside electron line is then 
inserted into the backward branch of the usual second 
order self-energy contribution. As we do not have an a 
priori knowledge about how to decorate the vertices with 
the Keldysh indices, we have checked all possibilities. It 
turns out that, at this order, there is only one single dia- 
gram that does not vanish for w > 0. It is shown in Fig. ||, 
a) and for small positive lu equal to n{Lo) — CUqV/lo, 
where C is a non-universal numerical constant (related 
to the suppression of the LDOS). For the FE problem we 
take F - Z?. 



a) 



b) 




FIG. 2. a) The only second-order diagram contributing 
to the TED above the Fermi energy. Solid lines correspond 
to electrons in the emitter (LL) and the dashed line to those 
on the right-hand-side of the contact. Crosses stand for the 
tunnelling vertices (7) and wiggly lines represent the interac- 
tions (f/o). b) Generic form of the diagram contributing to 
the TED for a; > 0. The four-particle vertex K, is denoted by 
a shaded square. 



In real systems 7 is supposed to be small but not 
Uf). Therefore our next objective is to develop a non- 
perturbative in Uq theory for the TED remaining at 
the 7^ order in the tunnelling. In the time domain the 
Green's function of interest g ^{ti, 12) is given by 

g~+{hM)^^{i>Hhml))=-^{Tcm^)^\tt)]) . 

where the superscripts indicate that the time variables 
ti and t2 lie on the time-ordered (T) and anti-time- 
ordered (T) parts of the Keldysh contour, respectively. 
A straightforward ^-matrix expansion in powers of 7 re- 
sults in (see Fig.||, b) 

Jc 

Disentangling the Keldysh indices one obtains four corre- 
lation functions that differ in the way the field-operators 
are combined under the r(T)-ordering operations. 



The best way forward is to derive a Lehmann- 
type spectral representation for the TED. When in- 
serting the complete set of states |A) at the break of 
the time-orderings and Fourier transforming back to 
the w-representation, we immediately find that of the 
four aforementioned correlation functions only one con- 
tributes to the TED for ui > 0: 



IC 



c 



IC = J2{mi'\h)4'\h)]\x){x\Tm^)^j{t^ 



Furthermore, inserting unity between the remaining field 
operators and computing the time-integrals we obtain 



(11) 



where = X),. a^^a^o[{Ef, - E^, + uj + iO) ^ + {Ei, + 
u! — iO)^^]. Greek indices count all excited states with 
energies E^^x,fj. and a^j, = {^liplv). Glearly all E^^'s are 
larger than the ground state energy Eq (which we have 
set to zero). Therefore we observe (though this fact is not 
relevant for the FE process) that the upper threshold for 
the TED, at the order 7^, is at w = (physically this 
is clear as the hot hole has only got so much energy to 
dissipate). We can push the spectral analysis further 
by noticing that if the T(T)-orderings in the /C-function 
were dropped then, due to analytic properties of time 
integrations, the overall result for n{LO > 0) would have 
been zero. Hence subtracting the un-ordered correlation 
function off and assuming the LDOS to be constant we 
can manipulate the time-integrals for the uj > TED 
into the expression: 

dt / dTi / dT2e''^(^i+^=) 

-00 t + ia Jo Jo 



X {{ij^{Tl), ij\0)}{^{t + Tl), ^{t + Tl + T2)}) 



(12) 



where {., .} stands for the anti-commutator and a ~ 1/Z? 
is the real-time cut-off. One advantage of this rep- 
resentation for the TED is that it is explicitly van- 
ishing in the non-interacting case as the field opera- 
tors then anti-commute at all times. The latter does 
not take place in interacting systems. For example in 
LLs, representation (^ leads to the following braiding 
relationship between the field operators, ip{ti)^{t2) = 
gi7rsign(ti-t2)/3^(^2)7/)(ii). Thus ouc possiblc interpreta- 
tion of the secondary TED is that it measures the degree 
of braiding induced by the interactions. 

So far our treatment has been general. Now we spe- 
cialise to LLs, in which case the four-point functions can 
be calculated explicitly. Using representation (|^) we ob- 
tain: 



JC[tiM-.hM) = (27r&) 2sgn(t3 - t2)sgn(t4 - ii) 
F{\t-i~t2\)F{\t^~t,\)F\Q) 



F{t2 - ti)F{t2 - ti)F{t3 - ti)F{t3 - ti) 



(13) 
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where F{t) = 1 ~ e"o(t+ia) ^nd eo = n/L. In the ther- 
modynamic hmit L — > oo we have to expand Eq. (^3|) in 
powers of eo and retain the leading term. The time in- 
tegrals involved can be computed in terms of generalised 
hypergeometric series but here we only give an in- 
tegral representation convenient to work with for small 
energies: 

fV-uj fji/g-i 
n{u;)^A{g,V) dE Fv{E + u), (14) 

where the spectral function is 

1*00 

Fv{p) = 2lme''''/3 / C^^^-'^e-^^^~P'>^ 

J a 

x^\l/g,0,-p( + ia), (15) 

and A{g, V) = T^l/g+l)^^/^ cos^7r/2g]/2nalr{l/g). 
Here 5* stands for the Tricomi confluent hypergeometric 
function. Expanding the latter on the branch cut fol- 
lowing [|l^ we obtain the asymptotic form of the TED 
in vicinity of the Fermi edge (at this final stage of the 
calculation we substitute the bias voltage by D): 

n{Lj)^CiiX + lfiu;/D)^, (16) 

where A — l/g — 2 and Ci is a non- universal constant, 
regular at A = — 1. In the limit of weak interactions, g 
1, the exponent A approaches the Fermi liquid result p7| . 
Eq. ( p^ ) can be regarded as consisting of two factors. The 
first factor is the universal 1 /to divergence inherent to all 
interacting systems. The second factor reflects power-law 
renormalisations normally occurring in the LLs and, in 
particular, contains the LDOS of the secondary electron. 
The latter object is suppressed for repulsive interactions 
{g < I). So also the TED singularity is suppressed. At 
gc = 1/2 the LDOS suppression effectively wins over and 
the singularity disappears. 

Because of the factorisation of the field operators, 
Eq. (H), the above results are valid for SWNTs as well. 
Using (0) we find that for SWNT Hi A = {3+l/K)/4:-2. 
Note that the critical value of the coupling now is Kc ~ 
1/5, which is actually within the experimental range (see 
above Eq. (^). Therefore wc do not make a specific pre- 
diction regarding the character of the singularity. Instead 
we think that both the divergent and the vanishing TEDs 
can be observed depending on the experimental setup. 

To conclude, in this paper we have put forward a the- 
ory for the FE effect in LLs and SWNTs. We have gener- 
alised Fowler-Nordheim relations and calculated the sec- 
ondary current that exhibits power-law singularities. We 
hope that this work will stimulate experiments on the 
FE from SWNTs in the relevant energy range, that is 
within a band ^ O.leV around Ep. More generally, we 
wish to stress that the secondary current is a fundamental 
effect due to a combined action of tunnelling and interac- 
tions. Measurements of TED spectra for different mate- 
rials could therefore reveal important information about 
the nature of electron correlations in particular media. 
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